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Listed  below  is  a  more  detailed  description  on  our  achievements. 

•  Dislocation  mobility  and  stability  in  nanocrystals  and  electronic  materials  are 
influenced  by  the  material  composition  and  interface  conditions.  Its  mobility  and 
stability  then  affect  the  mechanical  behaviors  of  the  composites.  In  the  JMPS  paper 
[1],  we  first  address,  in  detail,  the  problem  of  a  screw  dislocation  located  in  an 
annular  coating  layer  which  is  imperfectly  bonded  to  the  inner  circular  inhomogeneity 
and  to  the  outer  unbounded  matrix.  Both  the  inhomogeneity-coating  interface  and 
coating-matrix  interface  are  modeled  by  a  linear  spring  with  vanishing  thickness  to 
account  for  the  possible  damage  occurring  on  the  interface;  An  analytic  solution  in 
series  form  is  derived  by  means  of  complex  variable  method,  with  all  the  unknown 
constants  being  determined  explicitly.  The  solution  is  then  applied  to  the  study  of  the 
dislocation  mobility  and  stability  due  to  its  interaction  with  the  two  imperfect 
interfaces.  The  most  interesting  finding  is  that  when  the  middle  coating  layer  is  more 
compliant  than  both  the  inner  inhomogeneity  and  the  outer  unbounded  matrix  and 
when  the  interface  rigidity  parameters  for  the  two  imperfect  interfaces  are  greater 
than  certain  values,  one  stable  and  two  unstable  equilibrium  positions  can  exist  for  the 
dislocation.  Furthermore,  under  certain  conditions  an  equilibrium  position,  which  can 
be  either  stable  or  unstable  (i.e.,  a  saddle  point),  can  exist,  which  has  never  been 
observed  in  previous  studies.  Results  for  a  screw  dislocation  interacting  with  two 
parallel  straight  imperfect  interfaces  are  also  presented  as  the  limiting  case  where  the 
radius  of  the  inner  inhomogeneity  approaches  infinity  while  the  thickness  of  the 
coating  layer  is  fixed. 

•  In  the  Acta  Mechanica  paper  [2],  we  present  analytical  solutions  for  the  scattering 
of  antiplane  shear  wave  by  a  piezoelectric  circular  cylinder  with  an  imperfect 
interface.  We  first  consider  the  simple  case  in  which  the  imperfection  is  homogeneous 
along  the  interface.  Two  typical  imperfect  interfaces  are  addressed:  1)  mechanically 
compliant  and  dielectrically  weakly  conducting  interface;  2)  mechanically  compliant 
and  dielectrically  highly  conducting  interface.  The  expressions  for  the  directivity 
pattern  and  scattering  cross-section  of  the  scattered  antiplane  shears  wave  are  derived. 
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We  then  investigate  the  more  difficult  problem  in  which  the  imperfection  is 
circumferentially  inhomogeneous  along  the  interface.  A  concise  expression  for  an 
inhomogeneously  compliant  and  weakly  conducting  interface  is  derived  by  means  of 
matrix  notation.  Numerical  examples  are  presented  to  demonstrate  the  effect  of  the 
imperfection  and  the  circumferential  inhomogeneity  of  the  interface  on  the  directivity 
patterns  and  scattering  cross-sections  of  the  scattered  antiplane  shear  wave.  The 
circumferentially  inhomogeneous  interface  is  also  utilized  to  model  the  interface 
where  an  arbitrary  number  of  cracks  exist.  Results  show  that  when  every  part  of  the 
interface  is  rather  compliant,  large  low-frequency  peaks  of  the  scattered  cross- 
sections  as  shown  in  Fig.  1,  which  correspond  to  the  resonance  scattering,  can  be 
observed  no  matter  whether  the  interface  is  homogeneous  or  inhomogeneous.  The 
appearance  of  large  low-frequency  peaks  can  be  well  explained  by  estimating  the 
natural  frequency  of  the  corresponding  reduced  mass-spring  system  where  the 
cylinder  is  assumed  as  a  rigid  body.  Peaks  of  the  scattered  cross-sections  spanning 
from  low  frequencies  to  high  frequencies  can  be  observed  for  a  partially  debonded 


Fig.  1  Scattering  cross-sections  of  the  scattered  shear  wave  for  a  perfect  interface  (dashed 
lines)  and  for  a  rather  compliant  and  insulating  interface  a=0.ljU\/R,  /3=0  (solid  lines). 
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cylinder. 

•  In  the  JAM  paper  [3],  the  problem  of  a  functionally  graded  plane  with  a  circular 
inclusion  under  a  uniform  antiplane  eigenstrain  is  investigated,  where  the  shear 
modulus  varies  exponentially  along  the  x-direction.  By  introducing  a  new  function 
which  satisfies  the  Helmholtz  equation,  the  general  solution  to  the  original  problem  is 
derived  in  terms  of  series  expansion.  Numerical  results  are  then  presented  which 
demonstrate  clearly  that  for  a  functionally  graded  plane,  the  strain  and  stress  fields 
inside  the  circular  inclusion  under  uniform  antiplane  eigenstrains  are  intrinsically 
nonuniform.  This  phenomenon  differs  to  the  corresponding  homogeneous  material 
case  where  both  the  strain  and  stress  fields  are  uniform  inside  the  circular  inclusion. 

We  consider  an  infinite  FGM  in  the  x-y  plane,  and  assume  that  the  shear  modulus  /j. 
of  the  FGM  varies  exponentially  in  x-direction  as 

H  =  eipx/u o,  (1) 

where  /u0  is  the  homogeneous  shear  modulus  and  (5  the  gradient  factor  of  the  FGM. 

We  also  assume  that,  within  the  FGM,  there  is  a  circular  inclusion  r  -  -Jx2  +  y2  <  R 
which  undergoes  uniform  antiplane  eigenstrains  e\x  and  s\y .  The  boundary  condition 
along  the  inclusion-matrix  interface  r  =  R  is  assumed  to  be  fully  bonded,  and  can  be 
expressed  in  terms  of  the  out-of-plane  elastic  displacements  w(1)  inside  the  inclusion 
and  w(2)  outside,  as 

W(1)+w*=w(2), 

ew(1)  5w(2)  (r  =  *)  (2) 

dr  ~  dr  ’ 

where  w*  =  r(£^  -\s'zy)e'e  +r{s'’zx  +  \s]y)e~'e  is  the  additional  displacement 
corresponding  to  the  uniform  eigenstrains  s *a ,  s*  . 

We  now  introduce  a  new  function  <p  which  is  related  to  w  through  the  following 
relation 

w  =  e-px<p  (3) 
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It  is  easy  to  show  that,  in  terms  of  the  new  function^ ,  the  boundary  condition  (2)  can 
be  equivalently  expressed  as 


cp{1)  -<p°>  =  epxw* 


m  ( r  =  R) 

=  f3e^x  cos  0w  , 


9^<2)  8<f>  nBx 


(4) 


dr  dr 

where  <p(1,andy/2)are  within  and  outside  the  inclusion,  respectively.  They  satisfy  the 
following  Helmholtz  equations 


8V"  ,  8V" 


+2V_- av"=  o,  <r 


dxl  dy 


d2rnm  d2m{z>  i - 

M2)=o,  Jx^7>r 

dx  dy 


(5) 


In  view  of  Eq.  (5),  <p0)  and  <p{2)  can  be  expressed  in  terms  of  series  expansion  as 

•boo 

<pw  =  fJA(t)In{\P\rVne,  0<r  <  R  (6) 


<pm  =  2A»2)Kn( \PVVne,  r  >  R 


(7) 


where  I  and  K„  are  the  modified  nth-order  Bessel  functions  of  the  first  and  second 

n  n 

kinds,  respectively;  A ^  and  Ai2)  are  unknown  coefficients  to  be  determined. 


By  enforcing  the  boundary  condition  (4),  we  determine  the  unknown  expansion 
coefficients  in  Eqs.  (6)  and  (7)  as 


R 

'Ki\m 

-KM  R)~ 

U2)J 

ST 

-  c 

ST 

c 

1 

-  C 

* 

i 

J'MR) 

-i  MR). 

7„-i  W(4  -  ie*v )  +  /„+1  (/JR)(£l  +  'ys'zy) 

^|([/-2</h)+'.</»)](4  -i«;)+[/.(w+/„!o«)K4 +■«;>! 


(8) 


where  the  prime  (')  denotes  the  derivative  with  respect  to  the  variable  in  the 
parenthesis. 


•  In  the  IJES  paper  [4],  three-dimensional  Green’s  functions  are  derived  for  a  steady 
point  heat  source  in  a  functionally  graded  half-space,  in  which  the  thermal 
conductivity  varies  exponentially  along  an  arbitrary  direction.  We  first  introduce  an 
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auxiliary  function  which  satisfies  an  inhomogeneous  Helmholtz  equation.  Then  by 
virtue  of  the  image  method  which  was  first  proposed  by  Sommerfeld  for  the 
homogeneous  half-space  Green’s  function  of  a  steady  point  heat  source,  we  arrive  at 
an  explicit  expression  for  this  auxiliary  function.  With  this  auxiliary  function,  we 
finally  derive  the  three-dimensional  Green’s  functions  due  to  a  steady  point  heat 
source  in  a  functionally  graded  half-space.  Also  investigated  in  this  paper  are  the 
temperature  field  induced  by  a  point  heat  source  moving  at  a  constant  speed  in  a 
functionally  graded  infinite  space;  the  electric  potential  due  to  a  static  point  electric 
charge  in  a  dielectric  space  with  electric  field  gradient  effects;  and  two-dimensional 
time-harmonic  dynamic  Green’s  functions  for  homogeneous  and  functionally  graded 
materials  with  strain  gradient  effects. 


The  explicit  expression  of  the  three-dimensional  Green’s  functions  due  to  a  steady 
point  heat  source  in  a  functionally  graded  half-space  can  be  finally  derived  as 


The  temperature  field  induced  by  a  point  heat  source  moving  at  a  constant  speed  in  a 
functionally  graded  infinite  space  is 


(10) 

The  electric  potential  due  to  a  static  point  electric  charge  in  a  dielectric  space  with 
electric  field  gradient  effects  is 
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(11) 


</>= 


Q  l-exp(-/  ‘-\/x2  +  y2  +  z2 


An  e 


l 


x 2  +  y2  +z2 


It  is  of  interest  to  observe  that  Eq.  (1 1)  for  the  expression  of  the  electric  potential  due 
to  a  point  charge  in  a  dielectric  ceramic  with  electric  field  gradient  effects  is  still  finite 


(or  non-singular)  when  r  =  -Jx1  +  y1  +  z1  — >  0  (i.e.,  "when  approaching  the  location 
of  the  electric  charge).  More  specifically 

e 


lim^  = 
r~>°  An  el 

which  is  inversely  proportional  to  / .  Furthermore  it  can  be  proved  that 

e 


lim  *>«(-!)• 

<->o  dr” 


An  e{n  + 1)/" 


(/i  =  0,1, 2,3,...), 


(12) 


(13) 


where  =  <b  . 
dr° 


The  two-dimensional  time-harmonic  dynamic  Green’s  functions  for  homogeneous 
and  functionally  graded  materials  with  strain  gradient  effects  are 


w  =  ■ 


\n 


Hf{kxr)-Kfk,r) 


Inyfl^-Ak2!2 yi  L  2 
for  a  homogeneous  material,  and 

pi 2  (  f 

w  = - ,  — exp  ~—y 

lnl2y\  +  Ak2J2  V  2 


(14) 


1  jr 

-Hf(kir)-K0(k2r) 


(15) 


or 


w  =  ■ 


pi 


_exP 


2^/2Vl-4^2/ Vo 

for  a  functionally  graded  material. 


(16) 


•  In  the  AIAA  paper  [5],  we  address  in  detail  a  circular  inhomogeneity  with 
viscoelastic  interface  subjected  to  remote  uniform  antiplane  shear  stresses.  Both  the 
inhomogeneity  and  the  surrounding  matrix  are  assumed  to  be  elastic  and  quasi-static, 
while  the  interface  is  viscoelastic  modeled  by  a  linear  spring  and  dashpot.  Exact 
closed-form  solutions  for  both  the  Kelvin-and  Maxwell-type  viscoelastic  interfaces 
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are  obtained  by  means  of  the  complex  variable  method.  It  is  observed  that  when  the 
matrix  is  subjected  to  remote  uniform  shear  stresses,  the  stress  field  inside  the 
inhomogeneity,  although  time-dependent,  is  still  uniform.  The  derived  solutions  are 
then  used  to  predict  the  time-dependent  effective  shear  modulus  of  the  composite 
based  on  the  Mori-Tanaka  mean  field  approximation. 


By  using  the  Mori-Tanaka  mean-field  method,  the  time-dependent  effective  shear 
modulus  of  a  composite  containing  randomly  aligned  fibers  of  the  same  radius  with 
Kelvin-type  interface  on  the  xy  plane  can  be  derived  to  be 


M  =/V 


1  -c 

~Mi+fh  +  X<Jh-Mi)  2M  rrn| 

r  *+1,T 

.  (x+ViMi+Mi)  Cr  +  i Xm+^z) 

l  r  J. 

1  +  c 

+  2//, 

CAp 

.  Car  +  i Xth'+fh)  (z  +  iXh+Mi) 

l  r  J_ 

(17) 


where  //  is  the  time-dependent  effective  shear  modulus,  c  is  the  volume  fraction  of 
the  fiber.  It  is  seen  that  the  effect  shear  modulus  is  determined  by  the  shear  moduli  of 
the  fiber  and  the  matrix,  the  volume  fraction  c  of  the  fiber,  the  interface  rigidity  %,  the 
characteristic  time  y. 

The  time-dependent  effective  shear  modulus  of  a  composite  containing  randomly 
aligned  fibers  of  the  same  radius  with  the  Maxwell-type  interface  on  the  xy  plane  can 
be  derived  as 


1  -C 

1  exn! 

{  x1  T 

.  Or+iK/A  +  ^2) 

l  yOr  +  i)J. 

1  +  C 

",  ex_| 

{  x1  T 

.  Qr+1)Oi  +m2)  1 

l  r(z+i)J_ 

(18) 


It  is  added  that  the  method  presented  here  could  be  extended  to  address  the  interface 
described  by  a  standard  linear  solid  model,  which  combines  the  Maxwell-model  and  a 
linear  spring  in  parallel. 
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